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Self-induced transparency pulses propagating in a random medium embedded in a two-level system 
can transfer energy to localized Anderson states. This allows the onset of two-level laser-like action. 
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The fact that population inversion cannot be achieved 
at thermal equilibrium prevents the possibility of ob- 
serving laser action in a two level-system. [l| This sug- 
gests to resort to an out-of-equilibrium regime for a two- 
level laser, where light emission is amplified by an en- 
semble of atoms pumped at the resonant frequency; at 
variance with standard lasers, which adopt additional 
non-radiative decay levels that unavoidably and drasti- 
cally reduce the overall energetic efficiency (related is the 
case of Mollow gain [2]). A route is open by the well 
known Self-Induced Transparency process (SIT) 0-01 : a 
short optical pulse propagating in an otherwise absorb- 
ing medium such that the resonant two-level atoms are 
excited in its leading edge, and coherently decay in the 
trailing edge. This mechanism is described by a soliton 
solution of the relevant Maxwell-Bloch equations 1(3, II | , 
which, however, does not mediate a laser action, but al- 
lows the transport of a population inversion region in 
strongly absorbing media. Finding a way to localize such 
a region may lead to a two-level laser-like action. 
Here we propose to use a random system that sup - 
ports many overlapping electromagnetic resonances 12- 
18[, such that the probability of trapping the traveling 
population inversion is largely enhanced. The following 
issues are to be considered: (1) the effect of random- 
ness on the SIT soliton; (2) cavities are created by de- 
fects in an ordered structure with a band-gap; (3) the 
disorder-induced resonances must overlap with the SIT 
soliton spectrum. Points (2) and (3) lead to use a one- 
dimensional photonic crystal (PhC): a little disorder is 
indeed sufficient to create Anderson states in a known 
spectral region 19|, 2fj| |. 

For the point (1) above, we first consider the interaction 
of a SIT pulse with low-contrast refractive-index mod- 
ulations. We theoretically show that disorder doesn't 
destroy the pulse, which fluctuates around its mean po- 
sition with amplitude linearly growing with the propaga- 
tion distance. In a later part of the manuscript, through 
first-principle numerical simulations, we analyze a disor- 
dered PhC in the high-index contrast case; we numeri- 
cally find that at high scattering strength, the SIT pulse 
is trapped by coupling with band-gap Anderson states in 
proximity of the random structure surface 2l| . After the 
arrival of the SIT pulse, the random cavity starts emit- 



ting laser-like radiation with population inversion. 
We stress that non-solitonic pulses do not propagate in 
the absorbing medium containing the disordered PhC: 
using SIT is an unavoidable strategy to excite the dis- 
ordered system. In this respect, the interaction of SIT 
solitons with Anderson localizations provides a valuable 
route for a two-level laser-like action, while also being a 
fascinating regime for investigating the coupling between 
two forms of light localization respectively due to nonlin- 
earity and disorder 22- 2Sj|. 

Random SIT — We first address the role of a random 
perturbation 2(| 27| on a two-level 2k pulse by the re- 
duced Maxwell-Bloch (MB) equations. In a later section, 
we discuss the non-perturbative regime by resorting to 
first-principle numerical simulations by a parallel Finite 
Difference Time Domain (FDTD) code @. The MB 
equations with retarded time £ = t — z/v g and space z 
(v g is the velocity of the SIT pulse) read as 
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A = kB + S 1 (1) 



+ S 2 ,dtN = -{B*A + BA*) + S 3 , (2) 



being: A the (slowly-varying) electric field, B the atomic 
polarization and N the population invertion, 5*1.2,3 the 
perturbations, \x the detuning from resonant frequency 
wo, k = (4irujoNd 2 )/ch, d is the dipole moment, and c 
the vacuum light velocity. Let A = A s + A\ with A\ a 
perturbation and A s = (A s , B s , N s ) the soliton: 
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where w = ft (£ - X) and £ = exp{-ifi(£ - X) + i9}, 
fj, = rik/(/3 2 +fl 2 ). f3, X, 0, are the amplitude, position, 
detuning, and the phase of the soliton, being l/v g = 
l/c+ k/(2(3 2 + 2£1 2 ). The linearized system is 



(d,Xi) ■1 1 =£ 1 (A 1 ) + S P , 



(4) 
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where X\ — diag(l, 0, 0), S p is the perturbation and 
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Taking the SIT pulse resonant with the medium (fl = 
/i = 0), we derive equations for its parameters by first 
introducing the auxiliary functions [26| /x = 9x^4 s , /e = 
d g A s , fp = dpA s and f n = d n A s , and letting 



A 1 = f x 5X + fpSp + fg&e + faSQ + a R , 



(6) 



where <5A(z), S6(z), Sf3(z) and i5il(z) are the 
time-dependent perturbations to soliton parameters 
and clr is the radiation term (or = here- 
after, as it rapidly spreads and is absorbed. We 
let S p = {[V B {z)p + iV A {z)u]£, 0,0), such that 
(V A}B (z)V A , B (z')) = ( V a,b)$(z ~ where V a ,b are 
the electric field and polarization perturbation. Let- 
ting the adjoint functions [261] f s = ifp,fg = —ife, 

fa = -ifx, fx = ifn x , such that (f a ,f b ) = M a 5 a ,b, 
with a and b two parameters in (A, Q, 8, or j3); we have 
Af e =Afp = -l/(3/3 2 ) and Nx = Mj = -2tt/3 - 1/(3/3). 

Projecting over f k we get 



-k/p 2 sn+v A (z) (7) 



and (5/3 = S£l = 0, which gives 
(5X(zf 
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Eq. © describes the random fluctuations of SIT solitons 
in the limit of small uncorrelated disorder, and states 
that they grow linearly with the propagation distance, 
and decay when increasing of the soliton power: as the 
velocity of the soliton is reduced, its random walking be- 
comes more pronounced. This result describes an index 
perturbation with weak index-contrast, as shown below 
in figure [5^, figure [5^,b and figure for e r i and e r 2, 
in the comparison with FDTD simulations. Specifically, 
eq.© predicts that slower SIT solitons are those mostly 
affected by disorder, and hence they are expected to in- 
teract more effectively with Anderson states also in the 
high index contrast regime. As we numerically show in 
the following (figures [2}d, [3}; for e r 3, |4] and [5]), as the 
strength of disorder increases (when the theoretical anal- 
ysis reported above is not expected to the valid), this 
leads to localization and laser-like action. 
Numerical Results — Following Ref. (28j , we numerically 
solve the one-dimensional full MB equations by a paral- 
lel FDTD algorithm and study the evolution of the SIT 
pulse in a non-linear disordered SIT medium; we denote 
as P3 the population inversion of the FDTD equations 
(corresponding to N in Eq.© when adopting the rotat- 
ing wave approximation) and E is the real- valued electric 



field, whose slowly- varying complex envelope is A in (JTJ) - 
Note that we use a different notation for the field here 
because the FDTD equations are more general that the 
reduced MB equations © |2J} [28| . 

As sketched in figure [TJ the SIT pulse propagates into 
the grid from the left boundary (z = 0), from vacuum 
(z < 7.5/im) to a homogeneous SIT medium (atom den- 
sity N a toms = 10 24 m~ 3 ); its carrier frequency is resonant 
with the medium ojq = 27r/o, with fo = 2 x 10 14 s _1 , and 
the pulse duration t is chosen to satisfy the 2tt area theo- 
rem. The relaxation time of the density matrix equations 
are T 1= T 2 = 1.0 x 10" 10 s (T 1<2 » r). The initial pop- 
ulation inversion in the SIT medium is p^g = — 1, while 
in the dielectric layers p^o = and the two-level atoms 
are not present. We neglect absorption in the dielectric 
layers because much smaller than in the resonant two- 
level system. The input pulse, with peak Eq = (2h/d)/3, 
enters the two-level medium, initially set in the ground 
state P30 = —1, and generates the SIT soliton, which, 
after propagating in the homogeneous resonant region, 
interacts with the random structure as shown in Fig. [T] 
The total length of the structure is 200/im, the random 
active medium extends from 50/zm to 150/im, where we 
add a fixed degree of disorder by inserting slices of a non 
absorbing medium with relative permittivity e,.. The 
disordered structure is created by introducing random 
layers in the homogeneous two level medium. The de- 
gree of disorder is quantified by a parameter 7: the rel- 
ative permittivity distribution has a square-wave profile 
from 1 to e r (inset in figure la) given by the sign of the 
function sin(27rz/c? + 2^7^), with £ a uniform deviate 
in [0, 1] extracted in each point of the grid. If 7 = 
an ordered periodical structure is attained: the dielec- 
tric layers are equidistant, with constant width 200nm 
and period 400nm (d = 400nm) and displays a band-gap 
centered at ujq overlapped with the spectrum of the SIT 
soliton (see figure [TJd). For 7 = 0, increasing e r creates a 
one dimensional band-gap structure; for 7 > 0, increasing 
e r enforces the effect of disorder, and as shown in figure 
[2Jd in the absence of the two-level system; the localized 
states first appears in proximity of the band-edge of the 
ordered structure. We use this approach to selectivity 
create a distribution of Anderson states with a spectrum 
superimposed to that of the input SIT pulse. In the fol- 
lowing, the described disordered structure is embedded 
in a two-level system and SIT solitons are launched in the 
homogeneous region (z < 50/im) and then interact with 
the random system (z > 50/im). Non-solitonic pulses are 
rapidly absorbed and do not propagate. 

Figure [2] shows a temporal snapshot of the population 
inversion distribution in low index (e r = 1.5, 7 = 0) or- 
dered case and in the high index disordered case (e r = 11, 
7 = 1): in the former case, the wave propagates with 
limited distortion and the inversion population is trans- 
ported through the medium, in agreement with what ex- 
pected by the perturbation theory reported above (anal- 
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FIG. 1: (Color online) (a) The simulation region is 200/im 
long, the homogeneous two-level medium extends from 7.5pm 
to 50pm. The disordered two-level medium (e,- = 1) ex- 
tends from 50pm to 150/im; the initial population is equal 
to p3o = — 1 and randomly alternated by a dielectric non- 
resonant medium with a different refractive index (pzo = 0, 
e r = 11), as shown in the inset. A snapshot of the SIT pulse 
is shown before impinging on the structure; (b) SIT pulse 
spectrum for r = lOOfs (dashed line); spectrum of a trans- 
mitted short pulse without the resonant medium showing the 
photonic-band gap (continuous line, 7 = 0); transmitted spec- 
trum in absence of resonant medium with 7 = 0.5 showing the 
disorder induced resonances in the forbidden band. 
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FIG. 2: (Color online) SIT population inversion pz in (a) the 
low-index case (e r = 1.5, 7 = 0, analogous results for 7 > 0) 
and in (b) the high index case (e r = 11, 7 = 1), where surface 
Anderson states are excited. The insets show the correspond- 
ing E profiles. pz appears periodically modulated as pz = 
in the dielectric layers with e r > 1. 



ogous results for low index-contrast e r = 1 are obtained 
for 7 > 0, see figure [3Jd,c); in the latter case, the spatial 
distribution of the inversion population localizes in prox- 
imity of the input face of the sample. This shows that 



surface Anderson states are involved in the process [21|; 
as we will report in future work multiple SIT pulses can 
be employed to enlarge the trapped inversion region. 

Panels a,b in Fig. [3] show the trajectories and the corre- 
sponding standard deviation of the SIT pulse, calculated 



FIG. 3: (Color online) (a) SIT pulse trajectories calculated 
for 7 = 0.1, e r = 6, for three different initial pulse ampli- 
tudes £ i = 4.2186 x lOV/m, E 02 = 4-Eoi, E 03 = 8E 01 ; (b) 
corresponding standard deviations; (c) SIT pulse trajectories 
for three relative permittivities e r i = 2.5, e r 2 = 6, e r z = 11 
when Eo = -Em- Note that the curves for e r i,2 are cut (they 
proceed as straight lines) for comparison with the case e r z- 
The results are averaged over 100 realizations of the disorder. 



for different input peak values and for a fixed strength 
of disorder. In the low index contrast disordered case, 
following the previous theoretical analysis (see Eq. (jHJ), 
as the initial peak value increases, the soliton propa- 
gates faster in the structure, and the slope in Fig. [3^ 
increases and the corresponding fluctuations are reduced 
(Fig. [3)}) . Conversely, slow solitons perform a more pro- 
nounced random walk in the disordered structure. As the 
index contrast is increased (beyond the regime of validity 
of the perturbational approach reported above), this al- 
lows a localization processes as shown for e r ^ in Figure [3h 
displaying soliton trajectories for different strengths of 
the disorder (fixed 7 and increasing e r ). In the cases e r \ 
and e r2 the trajectories proceeds as straight lines, denot- 
ing a propagating pulse, conversely for e r s the trajectory 
bends and the soliton slows down and is trapped. In 
the other words, in the low-index contrast case, the soli- 
ton displays a weakly perturbed motion; as the strength 
of disorder increases (when increasing e r with fixed 7), 
its trajectory becomes more fluctuating, until the wave 
gets localized; correspondingly the system starts emit- 
ting laser-like radiation, as a population inversion region 
is formed in correspondence of an optical cavity. 

Figure 2] shows the output of the device for low and 
high index contrast for a fixed 7. For a small e r the 
input SIT pulse in is transmitted as shown in HJd, 
where the pulse exits at about 0.5ps. On the contrary, 
for high e r , the pulse is trapped in the structure and, af- 
ter a transient, the output (in Figf?];) corresponds to a 
laser-like emission from the the two-level system. Note 
the difference in the vertical scale: for Fig|4jD all the en- 
ergy is transmitted, while for FigHJ:, it is partially re- 
flected (see also figure E|) and partially trapped in the 
light-emitting localized states. In figure [S^,, we show the 
reflected temporal signal, comprising of the portion of 
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FIG. 4: (Color online) (a) Temporal profile of the input SIT- 
pulse with r = fOOfs and E = 4.2186 x 10 9 V/m; (b) trans- 
mitted pulse without disordered structure; (c) emitted signal 
for an highly scattering medium (e r = 11, 7 = 1) The hori- 
zontal scale in panels b,c is overlapped with panel a. 
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FIG. 5: (Color online) (a) Temporal profiles of the reflected 
electric field in the highly scattering case in Fi^iJ;, including 
the input pulse, its reflected fraction, and the subsequent two- 
level laser emission; (b) corresponding spectrum (continuous 
line) compared with the input SIT pulse (dashed line). 



the SIT pulse, which is not trapped in the disordered 
structure and is reflected (large peak around 0.3ps), and 
the subsequent emission from the population inversion 
trapped in the disorder medium; its spectrum in[5J} dis- 
plays characteristic peaks signaling the excited Anderson 
localizations with frequencies in the forbidden gap. The 
energy trapped in the localized modes decay with their 
characteristic lifetimes (of the order of 0.5ps). 
Conclusions — We investigated the interaction of a SIT 
soliton with Anderson localizations through theory and 
parallel Maxwell-Bloch simulations. We have shown that 
an increasing scattering strength of disorder progressively 
slows down the soliton, up to blocking the pulse in the 
random system. This process is accompanied by the ex- 
citation of modes in a disordered cavity in the presence of 
population inversion. This results into a two-level laser- 
like action. The interplay between various forms of light 
localizations, namely solitons and Anderson states, can 



hence lead to novel processes of light-matter interactions 
in complex systems. 
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